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Abstract 


A kinetics model for cellular repair and misrepair for describing mul- 
tiple radiation-induced lesions ( mutation-inactivation ) is coupled to a two- 
mutation model of initiation and promotion in tissue to provide a parametric 
description of tumor prevalence in the Harderian gland in a mouse . Dose- 
response curves are described for 7 -rays and relativistic ions. The effects of 
nuclear fragmentation are also considered for high-energy proton and alpha- 
particle exposures. The model described provides a parametric description of 
age-dependent cancer induction for a wide range of radiation fields. We also 
consider the two hypotheses that radiation acts either solely as an initiator 
or as both initiator and promoter and make model calculations for fraction- 
ation exposures from 7 -rays and relativistic Fe ions. For fractionated Fe ex- 
posures j an inverse dose-rate effect is provided by a promotion hypothesis 
using a mutation rate for promotion typical of single-gene mutations . 


Introduction 

An understanding of the deleterious biological 
effects of space radiation is needed before astro- 
nauts are subjected to prolonged exposures to the 
high charge and energy (HZE) ion component of 
the galactic cosmic rays (GCR). One of the primary 
concerns for deep-space flight is the expected can- 
cer risk from HZE ions. The nature of the ioniza- 
tions in tissue, including the track structure, from 
HZE particles is very different from any other type 
of radiation to which humans have been exposed, 
and thus the expected risk is largely unknown. Al- 
though experiments with HZE exposures with cell 
cultures are now quite numerous (Thacker, Stretch, 
and Stephens 1979; Yang et al. 1985; Kronenberg 
and Little 1989; Kranert, Schneider, and Kiefer 1990; 
Lett et al. 1989), only a few experimental stud- 
ies (Nelson et al. 1989; Ainsworth 1980; Burns and 
Albert 1980) of their mutagenic or tumorigenic po- 
tential in animal systems have been undertaken. 
The measurements of Fry et al. (1985) and Alpen 
et al. (1993 and 1994) for tumor prevalence in the 
Harderian gland of a mouse are the most useful of 
these studies because the effects on a single tumor 
type were considered, and dose-response curves for 
several ion species were studied using track-segment 
irradiations. 

A fairly common view in cancer research is that 
the transformation of a single cell will result in tu- 
mor formation (Pry and Storer 1987; Land, Parada, 
and Weinberg 1983; Renan 1990). The discovery 
of oncogenes and the mechanisms for their muta- 
tion has resulted in a widely held view of carcino- 
genesis for many tumor types as a multistep process 
(Renan 1990; Mitchel and Trivedi 1993), involving 


initiation, promotion, and progression. The initia- 
tion stage would include a set of mutations in the 
DNA of the cells produced by a carcinogen resulting 
in activation of one or more oncogenes (Renan 1990). 
The promotion stage describes the conversion of an 
initiated cell from a premalignant phenotype to a ma- 
lignant one, perhaps through the inactivation of a 
second type of gene called a tumor-suppressor gene, 
and, finally, progression is the stage of more aggres- 
sive tumor growth (Renan 1990). When radiation is 
acting as the carcinogen, the type of mutagens may 
vary widely. Photon irradiations are observed to in- 
volve point mutations (Renan 1990), and we should 
expect gross rearrangements and deletions in DNA 
to occur following HZE exposures. The GCR spec- 
trum imposes further difficulties because of the broad 
range of ion velocities and charge number that occur 
and because of the protracted exposure encountered 
on a long space mission. 

Several mathematical models of initiation and 
promotion have been developed for phenomenolog- 
ical descriptions of carcinogenesis, including mod- 
els of natural incidence of cancer (Moolgavkar and 
Knudson 1981) and radiation carcinogenesis follow- 
ing exposures to radon (Moolgavkar et al. 1990) or al- 
pha emitters (Marshall and Groer 1977). The works 
of Marshall and Groer (1977) and Moolgavkar et al. 
(1990) consider the effects of division and differentia- 
tion on initiated cells in order to provide a parametric 
description of the age-specific incidence. Models that 
undertake similar descriptions of cancer from GCR 
exposures must consider the continuum of carcino- 
gens entailed in the broad distribution of charge and 
energy of these ions and the success of HZE parti- 
cles for cell inactivation as observed in cell culture 


studies. Dose-response curves for tumor induction 
are often parameterized as (Upton 1986) 

/ = (c„ + c 1 D + c 2 £> 2 )e-(‘' D+w,I > (1) 

where I is the tumor occurrence, D is the absorbed 
dose, c 0 is the natural occurrence of cancer, b\ and b 2 
are the linear and quadratic coefficients, respectively, 
for high-dose saturation, and c\ and c 2 are the linear 
and quadratic coefficients, respectively. The expo- 
nential factor in equation (1) accounts for saturation 
or decrease in the occurrence at higher doses which 
is often attributed to cell inactivation. The useful- 
ness of equation (1) for HZE exposures is severely 
limited if methods are not available for determining 
the dependence of the parameter on radiation qual- 
ity. The track-structure model of Katz et al. (1971) 
has been successful in providing a parametric ap- 
proach for considering HZE effects with accurate pre- 
dictions made for an arbitrary ion species provided 
from fits to experimental data for 7-rays and a few 
ion types. Typically, the Katz model is considered 
only for acute exposures. More recently, a model of 
linear repair and misrepair kinetics has been devel- 
oped by Wilson, Cucinotta, and Shinn (1993) which 
includes the Katz action cross-section formalism and 
provides a description of temporal effects and the 
competition on a cell population between inactiva- 
tion and mutations. In this paper, we extend the 
kinetics model to include the assumption of initia- 
tion and promotion assumption in order to develop 
a parametric model of HZE carcinogenesis. 

The first purpose of this paper is to provide 
a model-dependent formalism that is the age- 
dependent analogy of equation (1) which should be 
useful for parametric descriptions of GCR effects 
based on the success of the Katz action cross-section 
model. This presupposes that the description of 
Katz for cell damage as described in the terminol- 
ogy of grain count, track width, and thindown car- 
ries over to cancer in animals from cellular stud- 


ies. A second purpose is to consider fractionated 
exposures, in particular, the role of cell inactivation 
and radiation induction of a second mutation of an 
initiation-promotion model. 

High-linear-energy-transfer (LET) ions have 
shown an inverse dose-rate effect in which fraction- 
ated or protracted exposures are often more se- 
vere than acute exposures (Ullrich, Jernigan, and 
Storer 1977; Ullrich 1984). The experiments of 
Ullrich suggest that for animal systems, protracted 
exposures with high-LET radiations participate in 
the promotion stage of carcinogenesis. Exposure 
with fission neutrons shows that the inverse dose- 
rate effect is tissue dependent as well as dose depen- 
dent (Ullrich, Jernigan, and Storer 1977). In cell- 
culture studies, synchronization experiments suggest 
that a sensitive window in the cell cycle exists (Hill 
et al. 1982; Brenner and Hall 1990). However, we 
do not know if the same explanation will be true in 
animal systems in which the length of the cell cycle 
is typically much longer than that in culture. Fi- 
nally, we consider model predictions for temporal ef- 
fects for the relative biological effectiveness of high- 
energy Fe and protons in which the effects of nuclear 
reactions or target fragments are included. In partic- 
ular, we consider any age dependence on the relative 
biological effectiveness. 

In the remainder of this paper, we first review the 
kinetics model of multiple-radiation-induced lesions 
and repair and misrepair of the lesions, including the 
introduction of the track-structure model. We next 
introduce a basic model of growth kinetics for ini- 
tiated cells and obtain expressions for tumor preva- 
lence. Dose fractionation is then considered in our 
model, including the possibility of radiation induc- 
tion of a mutation that we associate with the pro- 
motion stage of carcinogenesis. Finally, we discuss 
fits to the experiments for Harderian-gland-tumor 
prevalence and discuss model predictions for dose 
fractionation. 


Kinetics Model for Initiation 

\ kinetics model for cellular repair and misrepair has been developed by Wilson, Cucinotta, and Shinn (1993) 
that includes multiple-lesion formation and is based on first-order repair kinetics First-order kinetics models 
have been considered previously by Dertinger and Jung (1970) and Dienes (1966) for survival curves following 
photon exposures This development of Wilson et al. includes multiple-lesion types, such as those related to 
SSL mutation and utilizes the track-structure model of Kata et ah (1971) for mode, mg the 
lesion-formation rates for charged particles. We note that other forms of enzyme junet.es, 
order, second order, or mixed orders, can be considered. The use of first-order kinetics offers at leas 
simplicity of analytic solutions and provides a parametric framework. 
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The kinetics model assumes that nascent lesions are active chemical species produced in fast processes by 
radiation. The active chemical species, also denoted as substrates, are acted upon by enzymes in the cell 
and eventually repair to their original state or are misrepaired and left in a permanent damaged state. The 
fixation of the nascent lesions is assumed to occur over a time scale of minutes to hours and to follow first-order 
kinetics. As is well accepted in radiation biophysics, we have two distinct time scales, one for the initial events 
and another for the subsequential fixation of lesions. The uninjured population of cells at time t is denoted 
by no(0* The number of cells at time t wdth locus l damaged with a number i of lesions is denoted by nn(t). 
The production of these lesions by radiation is described by rate constants which will be dependent on 
radiation type (for example, the charge and velocity of an ion). The ability of the cell to repair damage leads 
to a rate of repair of the denoted a r / 2 , and if the active species are stabilized but left in a misrepaired state, 
a misrepair rate a m n occurs. (The rates for lesion formation and repair and misrepair are in units of inverse 
time.) The balance equations for the time development of cell populations in a single phase of the cell cycle 
are then given for the uninjured population (Wilson, Cucinotta, and Shinn 1993) as 

^0{t)=YL a rli n li( t )~ kn 0( t ') ( 2 ) 

li 

for a locus left in a misrepaired state as 

ni(t) =Yl a mli n li{ t ) + Y^ a rl'i' n l l 'i'{ t ) - kn l(t) ( 3 ) 

i Vi' 

and for the number of cells with l and i as 


i-l 

nu(t) = V] a rVi > n m ,(t ) + k H n 0 (t) + V n tj {t) - kn H {t) - a u n H {t) (4) 

Vi' 3 = 1 

In equation (4), n^y 2 <(t) denotes the cells with lesions at two loci which obey similar rate equations, and we 
define 

&li ~ a rli T Qmli (^) 

and 

k = ^2 k u ( 6 ) 

li 

We consider only the mutation or mutations at the loci associated with initiation with those cells left 
permanently fixed in this state denoted by n/(£). In order to proceed, loci associated with clonogenic death 
must be considered because the mutation phenotype must be expressed in a hereditary fashion. The radiation 
induction rates associated with clonogenic death are denoted by with 


k d = Y l k di ( 7 ) 

i 

The radiation induction rates associated w r ith the initiation mutation are denoted by kj t with 

kj = ^2 kj % ( 8 ) 

i 

The solution of rate equations (2)-(8) for acute exposures follows if we exploit the time scales for the fast 
radiation processes and subsequent fixation. Convergence in the dose range below 10 Gy is achieved for i < 3, 
and we assume a value of 3 for both clonogenic death and mutations. The solution for survival after repair is 
complete (Wilson, Cucinotta, and Shinn 1993) and is given by 


m d - 1 

n 0 (t)= n 0 (0)e~ kdU + ^ 


i'= 1 


^rdV_ 

a di' 



Q rli + OrdV 
a h + a di' 



( 9 ) 
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and for initiation by 


"/(<)= (e Mr - l) 


TTlfl 1 

e~ ktr n 0 (0) + £ ^n di ,(t r ) 
Q *' 


m d -l 


m d -l 

n Ii{tr ) + ^2 ~ E ^ Ln Iidi'( t r) 
i'=l ^ 


( 10 ) 


where m d = 3, mj = 3, t T is the duration time of the exposure, and where 


n 0 (t r ) = no(0) e ktr 
n dl (tr) = k d it r n 0 (0) e~ ktr 
n d 2 (tr) = (k d2 t r + ^k^tfjnoiO) e~ ktr 

ndi(tr) = (kd3U + ~k d2 k dl t 2 r + e~ ktr 

n Il(tr) = kjit r 7l 0 (0) e~ ktr 

nj 2 (t r ) = (fc/2<r + ^*n*r^«o(0)e -Wr 
2 

nildliU) = yk I1 k cn t 2 n 0 {0)e- ktr 

nildzVr) = {^k n k d2 t 2 r + —k^kr^fjnoib) e~ ktr 
n I 2 dl(tr) — (^^k d \kj 2 t 2 + -^k 2 ik d itf'jno(0) e ktr 


( 11 ) 


The possibility of radiation acting as both initiator and promoter suggests the treatment of two specific 
mutation types. In the experiments studying tumor formation in the Harderian gland of a mouse (Fry 
et ah 1985; Alpen et al. 1993 and 1994). animals are exposed at about 100 days and we should expect few 
cells already initiated for radiation to act upon. For fractionated exposures, the number of initiated cells will 
increase through radiation as well as division, and we will include a second mutation related to promotion 
which we describe below. We next discuss the treatment of track structure following Katz et al. (1971) in the 
kinetics model. 


Track-Structure Model for Lesion Formation 

The lesion-production coefficients in the kinetics model must include track-structure effects in order to 
describe HZE exposures. The model of Katz has been successful for many years in providing a parametric 
description of track structure and is used here to model the lesion-production rates. In the Katz model, 
biological damage from fast ions is assumed to be caused by secondary electrons (6-rays) produced along the 
path of the ion. The effects caused by energetic ions are correlated with those of 7-rays by assuming that the 
response in sensitive sites near the path of the ion is part of a larger system irradiated with 7-rays at the same 
dose level. The action cross section is the probability of single-particle (inactivation) activation or mutation 
and is calculated by integrating the 7-ray probability function over the radial path of the ion as 

0 ~ 2 tt / bdbP(b ) (12) 

JO 

where b is the radial distance from the ion to a sensitive site of characteristic size ao, T M is the maximum 6-ray 
range, and P(b) is a probability function for 7-ray response assumed to be of the multitarget or multihit form. 
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For example, in the multitarget model, 

P(6) = [l-e-^/ D o] m (13) 

where m is the target number, Do is the 7-ray radiosensitivity parameter, and D(b ) is the average dose at 
the sensitive site. The cross section calculated in equation (12) is observed to plateau at a value <tq, which 
is indicative of an effective damage area inside the nucleus. The cross section is observed to rise above the 
plateau value for stopping ions, which is referred to by Katz et al. (1971) as the track-width regime, and then 
to fall to zero, which is referred to as thindown. (Sec Katz, Dunn, and Sinclair 1985). For relativistic ions of 
moderate charge, a < <ro, and this is called the grain-count regime. In the track-structure model, a fraction 
1 — (cr/cro) of the fluence of the ion is assumed to be available to act through intertrack effects in a manner 
similar to 7-rays, and the gamma-kill dose of the ion is defined as 

M>-£) d (i4) 

with D 7 = 0 if a > 0*0? and D is the absorbed dose. In the grain-count regime, the action cross section is 
conveniently parameterized as 

a = a 0 (l- e- z ' 2 I^J (15) 


where Z* denotes the effective charge of the ion, /? denotes the velocity, and k is related to the parameters Dq 
and a q through 

K = P ° Q 0 X 1Ql1 (16) 

c 

where Do is given in Gy, ao is given in cm, and C is a constant that defines the average dose deposited in an 
extended target by an ion passing through that target. In the earlier work following the radial-dose model of 
Butts and Katz (1967), the constant C was set at 2. More recently, Chunxiang, Dunn, and Katz (1985) have 
considered a more accurate range-energy relationship for the maximum range of <5-rays. Using this model of 
radial dose, we find that C ~ 0.7, which effectively reduces the radius ao of the target by about 60 percent 
from earlier results using the Katz model. 

The lesion-induction coefficients of the kinetics model are matched to the Katz model by Wilson, Cucinotta, 
and Shinn (1993) through the choices 

k I2 t r = ojF (17) 

k d3 t r = o d F (18) 

where F is the fluence of the ion, 07 and 07 are the action cross sections for cell initiation and inactivation, 
respectively, and 


ki\t T = 6 


1/3 ^7/ 


k d \t r = 6 ' 


with all other values of k set to zero. The choices in equations (17)— (19) assume that 7-rays achieve only a 
single step, kji or k d \, whereas ions are capable of transition to the unrepairable states with a probability 
of kj 3 or fc d3 . For low-LET ions at high energies, the effects of target fragments must be included, and this is 
achieved by summing over the energy spectrum of all secondary ions as described by Cucinotta et al. (1991). 


Summary of Parameters for Initiation Kinetics 

• In order to clarify the number and meaning of parameters introduced in the kinetics model, we briefly 
summarize these. For describing the dose-response curves for 7-rays, a radiosensitivity parameter denoted 
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by Dq is used with a distinct parameter for initiation (Do/) and survival (D 0 <*)- Repair rates and efficiencies 
also occur in the kinetics equations. For acute exposures and for fractionated exposures with interfractionation 
times much longer than the time scale of repair (>1 day), only the repair efficiencies occur in the dose-response 
equations. In the model the cells that have sustained three or more nascent lesions are left fixed (repair 
efficiency is zero) as initiated cells or inactivated cells. This leaves tw'o repair efficiencies as parameters with 
values between zero and unity which are denoted by ct r d\/<*d\ and a rd2/ a d2 for cell survival and by a r /iA */ \ 
and ct r I2/ a I2 f° r initiation. The misrepair efficiencies are then determined as unity minus the repair efficiency. 
The repair efficiency for cells with one lesion is most important for fitting the dose-response curve, and this 
together with Dq determine the initial slope of the 7-rays (Wilson, Cucinotta, and Shinn 1993), which is zero 
for 100-percent repair efficiency. 

For describing the response to track-segment irradiations with charged ions, we also require action cross 
sections for both the cell inactivation and the initiation mutation. These are modeled by using the parametric 
track-structure model of Katz. Here, the response for any charged ion is determined from a knowledge of the 
radiosensitivity parameter for 7-rays for the identical end point, the average radial dose in a sensitive volume 
of radius ao, and an effective target area that encloses the sensitive volumes. By using equation (12), the 
action cross sections are then determined by fitting the values ao and <7 q to a data set as described by Katz 
et al. (1971). The cross sections for survival and initiation are distinct and result in two new parameters aQj, 
ao/ and <7o</, <tqi for each end point which are fit to dose-response curves. In practice, equation (15) is used 
for particles in the grain-count regime, with the result that the parameter k is used as the fitting parameter 
and ao is then determined by equation (16). For mixed-radiation fields (i.e., to include the effects of target 
fragments produced by high-energy ions), the contributions from all nuclear secondaries are summed to define 
an effective cross section as described by Cucinotta et al. (1991). 

Growth Kinetics and Tumor Prevalence 

In the two-mutation model of initiation promotion, the number of initiated cells must be specified as a 
function of tissue age. Here, the effects of division and differentiation of initiated cells are important for 
describing age-response curves (Moolgavkar and Knudson 1981). For describing the natural incidence, the 
kinetics equation for the normal cell population is assumed to be 

™o(0 = (-M/ ~ 00 + 70 )no(t) (20) 

where /x/ is the natural rate of initiation, /?q is the rate of cell loss, and 7 q is the rate of cell division with all 
rates in units of day' 1 . The initiated cell population (nj(t)) is determined by 

n(f)= fijn 0 (t)-{fx P + /?/ - Ji)nj{t) (21) 

where ftp is the natural rate of promotion, is the rate of cell loss for n/, and 7 / is the rate of division of 
initiated cells. The time rate of change of promoted cells is given by 

np(t)= up nj(t) (22) 

The rates for initiation, promotion, division, and cell loss in equations (20)-(22) may have some time 
dependence; however, they are assumed to be constant here. For the normal cell population, we will assume 
that most of the cells are quiescent (i.e., in the Go phase) and that losses are small, such that 

no(t) — se^^i 1 ^ s (23) 

where the initial population is denoted by no(0) = s. The solution for the initiated population with the initial 
condition nj( 0) = 0 is then 
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( 24 ) 


ni(t) = — [exp (7/ - fa - (ip)t - exp(-w«)] 

7/ - PI + /i/ - MP 

For 7/ — Pi fij or fip, the rate of growth is controlled by 7/ — /?/. 

The tumor prevalence is scored as the number of animals in which a neoplasm is found divided by the 
number of animals at risk at time t. In the kinetics model, the number of promoted cells is determined by 
equation (22) from which we define the hazard function (Marshall and Groer 1977) or rate of appearance of 
tumors as 

h(t)= h P {t - g) (25) 

where g is the minimum tumor growth time, which is some minimal time necessary for observing a tumor. The 
prevalence is then given by 

P(t) = 1 - exp (- / h{t) <ft) (26) 

We next include the effects of radiation induction of initiated cells for the case of acute exposures. The 
time scale of induction of lesions by radiation is certainly less than a fraction of a second. Enzymatic repair 
and misrepair of the lesions are observed to be complete in a few days, although for cancer induction, the state 
of knowledge is not well known. The kinetics of tissue growth will occur over a much longer time, perhaps 
many days, and we assume that the differences in time scales of these kinetic processes are such that they may 
be treated independently in a sequential manner. The effects of radiation and repair on no(t) and nj(t) are 
given by equations (8) and (9), respectively. By letting t r be the time of exposure, no(£ r ) and nj(t r ) be the 
number of normal and initiated cells immediately after exposure and repair are complete, respectively, and 
having irradiation occur early in the animal’s lifetime, we find that 

hp{t)= HJ£E£^h l (exp [(7/ - 01 - fipXt - tr} 

X exp(-fijtr) — exp(— /x/£)} + VP ™l(tr) exp|[7/ - /?/ - p P \t - t r J (27) 

The second term in equation (27) thus represents radiation-induced, initiated cells. An important question 
for modeling is whether radiation significantly modifies the growth of initiated cells. Here, we use the growth 
constant 7/ — /?/ from natural-incidence curves as a first estimate. Radiation will initially cause a blocking of 
progression through the cell cycle; however, the delay time should be small compared with the length of time 
elapsed before the observation of cancer, which is usually several hundred days. We note that experimental 
studies with 7-rays and neutrons show similar slopes for age versus incidence with the time of appearance 
shortened for high-LET radiations (NCRP 1990). We also note that the solutions given above for tumor 
prevalence versus time rely only on the combination 77 — /?/, and not on these parameters individually. 

Results for Fluence- Response Curves 

Tumor prevalence in the Harderian gland of B6CF 1 female mice after radiation exposure has been measured 
(Fry et al. 1985; Alpen et al. 1993) using pituitary hormones to increase the rate of expression. The exposures 
included 7-rays and several relativistic ions. Animals were exposed at an age near 100 days and were sacrificed 
at 600 days. The initial number of cells is about 5 x 10 6 per gland with a nuclear diameter of about 5.5 fim 
(according to a private communication with Fry in 1992), and we have assumed that about two-thirds of 
the cells are susceptible, with the result that the initial cell population is estimated at s = 2/3 x 5 x 10 6 . In 
figure 1 (see the dotted-line curve) we have fit the model of equation (26) to the data for natural incidence 
with pituitary isografts. The data arc from Fry et al. (1985) and from a private communication with Fry in 
1992, and the curves represent the choice for the parameters that are listed in table 1. The minimum growth 
time w r as set at 100 days with results not very sensitive to choices up to about 200 days. The limited amount 
of data was not sufficient to rigorously define the parameters; however, they are constrained to within about a 
factor of 2 for the model under study. 
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Table 1. Model Parameters for Harderian Gland Tumors 


(a) Natural-incidence parameters 


/i/, per day lx 10 -7 

/ip, per day lx 10~ 7 

7/ — j3[, per day 6 x 10~ 3 


(b) Radiation-induction parameters 


End points 

CTO* cm 2 

K 

m 

A)> Gy 

Survival 

3.2 X lO -7 

550 

3 

3.2 

Initiation 

7.6 x lO -10 

480 

3 

148.0 


(c) Repair efficiencies 



Values of a n /a, for 

- 

End points 

i = 1 

i = 2 

Inactivation 

0.999 


0.5 

Initiation 

0.995 


0.5 


Model 


600-MeV/amu Fe at 64 cGy 

600-MeV/amu Fe at 10 cGy 

Natural incidence 



Age, days 

Figure 1. Prevalence plotted against age of mice at risk for 
Harderian glad tumors using pituitary' isografts. 

In figure 2, fits to the data of Alpen et al. 
(1993) for tumor prevalence using equations (25)— 
(27) and (10) are shown versus particle fluence. An 
LET value of 0.23 keV/fim is assumed for 7-rays. A 
summary of the ion types and their energies and lin- 


ear energy transfers are given in table 2. Calculated 
values are also given in table 2 using equations (12) 
or (15) for the action cross sections for cell inacti- 
vation and initiation with the fitted parameters for 
radiosensitivity and action cross sections in table 1. 
The repair efficiencies are also listed in table 1. The 
error bars in figure 2 denote the standard deviations 
of the prevalence reported by Alpen et al. (1993). 
Overall, the agreement with the data is good for all 
but two data points fitted to within the experimental 
95-percent confidence intervals (Alpen et al. 1993). 

For the *H and 4 He exposures, the effects of target 
fragmentation were included following the method 
described by Cucinotta et al. (1991), and their contri- 
butions to a are listed in parentheses in table 2. The 
target fragments represent a substantial increase in 
the prevalence as compared with the 7-ray response 
for these low-LET ions, as can be seen in figure 2(a). 
The 4 He response, however, is underpredicted at the 
higher fluences. The present model offers no expla- 
nation for the differences seen between 1 H and 4 He at 
large fluences. The model predicts and the data sug- 
gest a turn-down in the prevalence at large fluence 
for 4 He which we attributed solely to cell killing, 
thus reducing the number of target cells available for 
initiation. 

For the Nb exposure, the cross section is calcu- 
lated through integration over the radial distance 
using equation (12) because the fit equation is not 
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Fluence, Jim -2 


(a) Low-to-medium LET radiations. 



(b) High-LET radiation. 


Figure 2. Prevalence of Harderian gland tumors at age of 600 days plotted against particle fluence. Experimental data are taken 
from Alpen et al. (1993). Error bands represent standard deviations. 
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Table 2. Radiation Types and Action Cross Sections of Tumor Prevalence in Harderian Gland 


Radiation 

type 

Energy, 

MeV/amu 

LET, 
ke V / /im 

07 , cm 2 
(a) 

cm 2 

(a) 


250 

0.4 

1.3 x lO " 16 (5.0 X 10~ 14 ) 

3.6 x 10 -14 (1.9 x 10 -11 ) 

4 He 

228 

1.6 

9.7 x lO " 15 (1.0 x 10" 13 ) 

2.7 x 10 - 12 (3.8 x 10- n ) 

20 Ne 

670 

25 

1.5 x 10 -11 

4.5 x 10 -9 

56 Fe 

600 

193 

5.4 x lO -10 

2.0 x 10 -7 

56 Fe 

350 

253 

6.5 x HT 10 

2.5 x 10 “ 7 

93 Nb 

600 

464 

1.0 x 10~ 9 

5.0 x 10 “ 7 


“Values in parentheses represent target- fragment contributions. 


Table 3. Cellular-Response Parameters for Survival in Mammalian Cell Lines 


Cell type 

(Tq, cm 2 

K 

m 

Do, Gy 

Harderian gland 

3.2 X 10 “ 7 

550 

3 

3.2 

C3H10T1/2 

5.0 x 10 “ 7 

750 

3 

2.8 

V-79 

4.28 x 10 ~ 7 

1100 

3 

1.82 

T-l kidney cells 

6.7 x 10 -7 

750 

2.5 

1.7 

HeLa 

5.6 x 10 - 7 

1100 

3 

3.7 

Mouse bone marrow .... 

4.2 x lO " 7 

500 

2.5 

0.9 


accurate in the track-width regime. We note that 
no reduction in carcinogenic potential is expected 
from the highest LET ions for the ions under study 
in the present model because of their relatively high 
velocities. For stopping ions which have similar LET 
as the Nb beam, the model would predict such a 
reduction. The parameters obtained for the initiation 
cross section (table 2) estimate an effective area <tq 
slightly larger than that seen in mutation studies 
(Tsuboi, Yang, and Chen 1992) or transformation 
studies in vitro (Yang et al. 1985) which suggests 
that several genes are able to act as initiators. 

Transplantation studies of Harderian gland cells 
from CBA/Cne mice into the fat pads of isogenic re- 
cipients were studied by Di Majo et al. (1986) where 
in vivo survival curves were measured following X-ray 
irradiation. The age of the mice at exposure was ap- 
proximately the same as the ages used by Fry et al. 
(1985) and Alpen et ah (1993). Although the mice 
are of a different strain and the experiments of Alpen 
et al. (1993) used a 7 -ray source for the reference ra- 
diation, we compare their survival measurements in 
figure 3 with the model survival curves that result 
from our model as given by equation (9). The agree- 
ment is good, but the relevance is uncertain for the 
reasons stated. Although in vivo measurements (us- 



Figure 3. Survival plotted against absorbed dose in Harderian 
gland tissue. Experimental data are taken from Di Majo 
et al. (1986) for CBA/Cne mice using X-ray irradiations. 
Error bands represent standard deviations. 

ing charged particles) for this tissue are not available, 
we note that the cellular-response parameters for sur- 
vival listed in table 1 are similar to those fitted by the 
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Katz model for many mammalian cell lines, which, 
for comparison, are listed in table 3. Cell killing has 
an important effect on the tumor prevalence mea- 
sured at intermediate and large doses for heavy ions 
and cannot be ignored in model predictions of tumor 
induction. 

Age Effects and Role of Pituitary 
Isografts 

In figure 1, we have calculated the tumor preva- 
lence as a function of age for the 600-MeV Fe ex- 
posure at doses of 10 and 64 cGy. Also plotted are 
data for a fission neutron (fn) exposure with a mean 
neutron energy of 0.85 MeV at a dose of 64 cGy. 
Calculations were not performed for the fn exposure 
at this time because of the detailed transport analy- 
ses required. The similarity in the time development 
for Fe and fn is quite noticeable, with the slope of 
the model calculations for the prevalence curve for Fe 
being more rapid than that of the experiments with 
fission neutrons. We noted that an identical slope 
could be achieved by using a decrease in the growth 
rate with a corresponding increase in the spontaneous 
promotion rate. The comparison in figure 1 sug- 
gests that using the value 7/ — /?/ = 6xl0~ 3 per day 
fitted to the natural-prevalence curve for radiation- 
induced tumors is fairly accurate for the system un- 
der study. The use of pituitary hormones is expected 
to increase the growth parameter 7 j — (3j over the 
natural rate in the two- mutation model. Exposures 
were also performed for Fe (Alpen et al. 1993) and fn 
(Fry 1981) without the use of pituitary isografts. In 
both cases, a decrease appeared in occurrence; how- 
ever, for the Fe exposures, the decrease was small. 
We fit our model to the data of Grahn, Lombard, and 
Cranes (1992) for the natural prevalence of Harder ian 
gland tumors without isografts and found a value of 
7/ - ^ 2.7 x 10 -3 per day if we keep pj and \ip 

fixed as in table 1. The data for Fe without isografts 
can be fit in our model with 7 /—/?/ = 3.5 x 10~ 3 per 
day if all other parameters are unchanged. 

In figure 4, we show calculations of the relative 
biological effectiveness (RBE) for at 250 MeV 
and 5 °Fe at 600 MeV/amu as a function of age for an 
excess prevalence of 3 percent. The calculations sug- 
gest that for acute exposures, the RBE is sensitive 
to age with an RBE increase with age due to the re- 
duced role of cell inactivation in achieving a 3-percent 
excess prevalence at the later ages for these ions. 

Dose Fractionation and Promotion By 
Radiation 

"Many experimental studies in animals as well as in 
cell culture observe an enhancement in oncogenic re- 



Age, days 


Figure 4. RBE at 3-percent excess prevalence for Harder ian 
gland tumors plotted against age of mice exposed at 
100 days. Model calculations are for 56 Fe and ] H. 


sponse for protracted or fractionated exposures with 
high-LET radiations. This is in contrast to photons 
where a sparring effect is observed as the norm. Sev- 
eral of the possible explanations for the enhancement 
or inverse dose-rate effect include a sensitive phase in 
cell cycle, the effects of reduced cell killing for pro- 
tracted exposures els compared with acute exposures, 
the role of repopulation for tissue systems, and, fi- 
nally, the possibility that radiation will act as a pro- 
moter of initiated cells. The recent experiments of 
Miller et al. (1990) using synchronized C3H10T1/2 
cells establish that a sensitive phase exists for the 
transformation of cell cultures. For in vivo carcino- 
genesis, a large amount of data (Upton 1986; Ullrich, 
Jernigan, and Storer 1977; Ullrich 1984) performed 
under varying exposure conditions suggest that sev- 
eral factors contribute to the inverse dose-rate effect. 


In principle, the kinetic equations in equa- 
tions (2)— (1 1) could be extended to include the cell 
cycle (Wilson, Cucinotta, and Shinn 1993), as well 
as a second mutation type corresponding to promo- 
tion which could be interpreted as the activation of 
a tumor-suppressor gene. Such an approach would 
be cumbersome because of the need for a numerical 
solution of the resulting differential equations and be- 
cause of the large number of parameters that would 
be required. Instead, we will consider fractionated 
exposures with interfractionation times longer than 
a few days in which analytic solutions are possible 
without the addition of many new parameters. This 
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allows us to consider several of the proposed expla- 
nations of the inverse dose-rate effect in the present 
model. 

The roles of cell killing and repopulation are 
considered for a fractionated exposure separated by 
time intervals of a few days or more by sequential 
solution to equations (2)-(4) and (20)-(22). Here, 
the number of initiated cells from radiation at time t 
after N fractions for the case of no radiation-induced 
promotion is 

N 

n l( t )='^2 n d t i) ex p[(7/ -01- - UVi ( 28 ) 

i= 1 

where nj(ti) is the number of cells initiated by ra- 
diation in the zth fraction and fj is the fraction of 
cells remaining after the zth fraction. To consider 
an upper bound on fractionation effects from possi- 
ble repopulation between exposures, we also consider 
setting fi equal to unity which corresponds to the 
case of full repopulation. 

We will also consider a second mutation event 
corresponding to the promotion of initiated cells in 
fractionated exposures. This type of assumption be- 
comes important for older ages in which the num- 
ber of cells initiated spontaneously will be relatively 
large. Here, the kinetic equation solution follows 
from equations (2)- (4) with the introduction of a 
second mutation type that is active only in the ini- 
tiated cells. If we assume that the populations nj 
and np are always small compared with tt-q, we find 
that the number of radiation-induced promoted cells 
in an TV-fraction experiment is given by 

( 29 ) 

where np{ti)/nj(t) represents the fraction of initi- 
ated cells promoted by the zth exposure and is given 
by equation (10), with the lesion formation rates and 
repair rates for promotion used instead of those for 
initiation. We estimate these rates in the following 
comparisons. 

Figure 5 shows four calculations for fractionated 
exposures separated by 1 week of 600-MeV / amu Fe 
and y-rays versus absorbed dose. In figure 5(a), we 
use equation (28) with no allowance for repopulation 


between fractions. Two effects are observed: (1) the 
7-rays show considerable sparring with increasing 
number of fractions, and (2) for Fe, a 24-week frac- 
tionated exposure also shows considerable sparring 
that is due to the initiated cells in the later fractions 
having insufficient time to divide before the sacri- 
fice at 600 days. In figure 5(b), we allow for full 
repopulation between exposures. Again, the 7-rays 
show sparring. For Fe, an inverse dose-rate effect oc- 
curs for the two- and six-fraction schedules; however, 
for the 24-week schedule, the increase is seen only 
above 0.8 Gy w r ith sparring again at the lower dose 
levels. The inverse dose-rate effect seen in fission neu- 
tron exposures for many tumor types shows a con- 
tinued increase for even 60- week schedules (Grahn, 
Lombard, and Cranes 1992). Because we also do 
not expect full repopulation to occur, we expect that 
figure 5(b) does not account for the enhancement 
anticipated. 

In figure 5(c), we include a second mutation rate 
from equation (29) and have all lesion formation 
parameters set equal to those used for initiation. 
A large enhancement is now seen for the fraction- 
ated exposures with a decrease for 7-rays again. 
With the fractionation combined with the mutation 
rate, we have assumed that the second mutation 
accounts for an inverse dose-rate effect. As a fi- 
nal estimate, we have considered calculations in fig- 
ure 5(d) in which the promotion lesion is assumed 
to be similar to the HGPRT (hypoxanthine guanine 
phosphoribosyl transferase) mutation as estimated 
in our mutation model from the data of Tsuboi, 
Yang, and Chin (1992) and Thacker, Stretch, and 
Stephens (1979). Here a$ p = 0.9 x 10~ 10 cm 2 and 
Dq p « 1000 Gy. For this very specific gene muta- 
tion, the enhancement is seen, although it is less sig- 
nificant than it would be if the promotion mutation 
rate were near that of initiation . 

In space, the dominant radiation component is 
high-energy protons in which a flux of about 2 x 10 s 
protons/cm 2 /year is expected. Because these pro- 
tons will produce an appreciable number of high-LET 
secondaries, we consider if an inverse dose-rate ef- 
fect will occur for protracted exposure. In figure 6, 
we show calculations for 250- MeV protons using the 
model of figure 5(c). A large enhancement is seen 
above doses of 0.5 Gy, which is close to the expected 
dose on an extended space mission. 
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Figure 5. Model calculations of dose fractionation for Fe at 600 MeV/amu and 7 -rays. Interfractionation times arc 1 week. 



Figure 6 . Model calculations of dose fractionation similar to 
figure 5(c) for 250-MeV protons. 


Concluding Remarks 

By using the hypothesis that carcinogenesis in 
mice occurs through two mutational steps, we have 
developed a parametric model of radiation carcino- 
genesis for charged particles. The number or type of 
mutations required for cancer induction is not well 
known and is certainly not unique. By assuming 
a two-mutation model with clonal expansion of ini- 
tatied cells, the age dependence of natural occurring 
tumors can be fitted, and the possibility of radiation- 
induced promotion can be explored. The model pro- 
ceeds from the kinetics of lesion formation and repair 
and misrepair for the mutation and survival of cells. 
The use of linear-repair kinetics provides an analytic 
framework to consider dose-rate effects. However, 
many important questions regarding the kinetics of 
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enzymatic repair have not been considered, and these 
may become important in extrapolating a parametric 
model to low fluences. Track- structure effects have 
been introduced into our model through the use of 
the radial-dose formalism of action cross sections de- 
veloped by Katz. The resulting parameterizations of 
action cross sections for mutations and inactivation 
as a function of charge and velocity of an ion allows 
for predictions for any monoenergetic or mixed field 
of radiation for which the particle- fluence spectrum 
is known. 

A two-mutation model of the natural incidence 
of carcinogenesis requires rates for spontaneous pro- 
duction of the first and second mutations, as well 
as the rate of clonal expansion of the initiated cells 
that carry the first mutation. We have estimated 
these rates by using the natural-incidence curves for 
Harderian gland tumors in mice that have received 
pituitary isografts. The kinetics of radiation-induced 
mutation were coupled to the model of the natural 
incidence of cancer. For mice having acute exposures 
early in their mature life, we have assumed that the 
rate of expansion of initiated cells is close to the spon- 
taneous rate. This assumes that the rate does not 
change appreciably with radiation type or damage 
level and that radiation-induced blocking of the cell 
cycle has only a small effect on the expansion sev- 
eral hundred days after the administration of radia- 
tion. The resulting model was fitted to dose-response 
curves for Harderian-gland-tumor prevalence in mice 
near 600 days in age. The cross sections for inactiva- 
tion and the mutation associated with the initiation 
event determined from our fits are of the same order 
of magnitude as those observed in many experiments 
with cell culture. For high-energy protons and alpha 
particles, the addition of the effects of the target frag- 
ments produced in nuclear reaction accounted for the 
increase in tumorigenic potential seen at low dose as 
compared with 7 -rays. The large differences seen in 
the experiments between protons and alpha particles 
at high dose could not be explained in our model. 

An enhancement in onconogenic effect following 
protracted exposures to high-linear-energy-transfer 
(LET) radiation has been observed in many stud- 
ies in animals and cell culture. By using our ap- 
proach, we have considered several factors that could 
lead to such an effect, including cell killing, repopu- 
lation, and radiation-induced promotion of initiated 
cells. The effects of repopulation were seen to lead 
to an enhancement for a small number of fractions 
for relativistic iron nuclei; however, for a large num- 
ber of fractions, the enhancement was not seen be- 
cause of insufficient time for expansion of initiated 
cells. The addition of a second mutation induced by 


radiation associated with the promotion of initiated 
cells also leads to an inverse dose-rate effect in the 
present model. The enhancement using this assump- 
tion would be very large if the action cross section 
for the second mutation was about the same as that 
of the first mutation, or it would be more modest 
if the cross section was close to the observed muta- 
tion rates in mammalian cells for the HGPRT locus. 

In all cases considered, the effects of 7 -rays are re- 
duced through dose fractionation. In contrast, an in- 
verse dose-rate effect is seen for high-energy protons 
when radiation is assumed to act as a promoter be- 
cause of the high-LET component of their effect from 
nuclear reactions. This effect could have important 
consequences for space radiation protection. 
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